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We propose a novel scheme for high fidelity photonic controlled-phase gates using Rydberg blockade in an 
ensemble of atoms in an optical cavity. The gate operation is obtained by first storing a photonic pulse in the 
ensemble and then scattering a second pulse from the cavity, resulting in a phase change depending on whether 
the first pulse contained a single photon. We show that the combination of Rydberg blockade and optical 
cavities effectively enhances the optical non-linearity created by the strong Rydberg interaction and makes the 
gate operation more robust. The resulting gate can be implemented with cavities of moderate finesse allowing 
for highly efficient processing of quantum information encoded in photons. As an illustration, we show how the 
gate can be employed to increase the communication rate of quantum repeaters based on atomic ensembles. 
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Large bandwidth, fast propagation and the non-interacting 
nature of photons, make them ideal for communicating quan¬ 
tum information over long distances H]. In contrast, strong 
photon-photon interactions are desirable for processing of 
quantum information encoded in the photons, especially if 
both high fidelity and high efficiency are needed. To satisfy 
these requirements one needs a highly non-linear medium. 
Typically, the strength of photon-photon interactions medi¬ 
ated by a non-linear medium is very weak at the single-photon 
level where photonic quantum logic gates are operating El. 
As a consequence, the implementation of photonic quantum 
gates remains an unsolved challenge and requires novel means 
of efficient light-matter interaction. To enhance light-matter 
interactions, a viable solution is to use ensembles of atoms, 
e.g., configured for electromagnetically induced transparency 
(EIT) 0; this can be further improved by placing the ensem¬ 
ble in an optical cavity, but these ensemble based approaches 
do not increase the essential photonic non-linearity. In recent 
years, there has been intense efforts to realize light-matter in¬ 
teractions via, non-linear interactions in a variety of medium, 
ranging from atoms BGl and atom like systems ifTTmrl to 
superconducting qubits HMnl. 

A promising approach towards creating strong quantum 
nonlinearities is to exploit excitation blockade in Rydberg EIT 
systems iHMa. Several quantum effects like strong optical 
non-linearities and control of light by light ll22l - E8l . determin¬ 
istic single-photon sources 1^ . and the generation of entan¬ 
glement and atomic quantum gates If30l434ll have been inves¬ 
tigated. The strong nonlinearity originates from the fact that 
the Rydberg interaction prevents multiple excitations within 
a blockaded radius rf, ||35]|36l. The ensemble then behaves 
as a two-level superatom consisting of Nb atoms within a ra¬ 
dius Tb llTSlI^ . If the optical depth db corresponding to the 
superatom is sufficiently large, db ^ 1 ll22ll . a strong optical 
nonlinearity at the single-photon level can be achieved in the 
EIT configuration lfT8l - l25]l . Reaching such an optical depth is. 


however, challenging, which limits the effectiveness of pho¬ 
tonic quantum gates. 

An enhanced optical nonlinearity was recently demon¬ 
strated by placing the ensemble inside an optical cavity 1261, 
but a direct application of this nonlinearity for quantum gates 
is non-trivial since the outgoing optical modes are highly dis¬ 
torted and entangled by the interaction llJ7l439l . In this let¬ 
ter, we propose a novel scheme for achieving a high fidelity 
photonic controlled-phase (CP) gate with a Rydberg EIT en¬ 
semble trapped inside an optical cavity of moderate finesse. 
In our scheme, the photons are incident at different times thus 
avoiding the problem of mode distortion while still allowing 
the cavity enhancement of the interaction. The use of a cav¬ 
ity has several major advantages compared to ensembles in 
free space, since it enhanced light-atom coupling in the en¬ 
semble and also effectively increases the non-linearity. In our 
proposal, the parameter characterizing the Rydberg blockade 
is Cb ~ d^db, where Cb — NbC, with C ^ 1 being the sin¬ 
gle atom cooperativity, and F is the cavity finesse. Hence the 
effect of the Rydberg interaction is increased by the cavity fi¬ 
nesse F, whereas the low value of C is compensated by a high 
value of Nb- In addition, the cavity is also useful for con¬ 
trolling the mode structure thereby enabling high input-output 
efficiencies iQi. We show that the proposed gate can have 
a promising (heralded) error scaling as l/C^, and demon¬ 
strate how it can be used to improve quantum repeaters based 
on atomic ensembles even for moderate interactions strengths 
Cb 10. The proposed CP gate can thus be directly inte¬ 
grated into quantum communication circuitry thereby provid¬ 
ing a building block for future quantum networks. The Ryd¬ 
berg interaction ll^ for our proposal can either be long range 
dipolar or van der Waals interactions, but for simplicity, we 
only consider the latter. 

We first outline the basic idea of our gate, which goes along 
the line of Ref. BTl . except that the single trapped atom is 
replaced by a Rydberg ensemble. In contrast to Ref. HD, 
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and many others, we thus do not require the strong-coupling 
regime of cavity QED, and can work with cavities of moderate 
hnesse. This enables input-output efficiency near unity since 
the cavity losses can be completely negligible compared to the 
mirror’s transmission M- For simplicity, we first describe 
the operation for single-rail qubits where a qubit is encoded in 
a photon pulse containing a superposition of vacuum 10) and 
a single photon |1). Later, we generalize it to a more useful 
dual-rail encoding where the qubit is encoded as a photon in 
one of two possible modes. 

In the single-rail version outlined in Figs. [TJa) and (b), a 
first photon pulse is stored in a cavity containing a Rydberg 
FIT ensemble ll42l . Here a classical driving field from an ex¬ 
cited state |e) to a Rydberg state |r) enables the storage of 
incoming photons in |r) through the interaction of the cavity 
held with the transition from the ground state \g) to |e). The 
excitation in state |r) is then transferred to another Rydberg 
state |r') by a microwave pulse so that the ensemble contains 
a single atom in state |r') if the hrst pulse contained a single 
incoming photon. The second pulse is then incident on the 
cavity. If the hrst pulse contained vacuum |0), the second 
pulse is scattered under Rydberg FIT conditions and leaves 
the cavity with the same phase. If the hrst pulse contained a 
photon, the atom in \r') shifts the position of the state |r) in 
the remaining atoms. As we will show, this prevents the sec¬ 
ond pulse from entering the cavity resulting in a phase hip on 
the 11) component of the second pulse. This evolution thus re¬ 
alizes a CP gate which, together with single qubit operations, 
is universal for quantum information processing. 

We now present a theoretical treatment to evaluate the per¬ 
formance of the CP gate. The initial state of the single photon 
pulse can be expressed as / c?a;(/)(a;)aj,e““*|0), where 
is the normalized pulse shape, is the one-dimensional held 
operators satisfying the standard bosonic commutation rela¬ 
tions and |0) denotes the vacuum of all the optical modes. 
The frequency integrand w of the incoming photon is refer¬ 
enced to the cavity frequency Wc, which in turn is nearly res¬ 
onant to the |e) —^ \g) transition (see Fig. [^a)). The cavity 
is assumed to be one-sided with a standing-wave held. The 
dynamics of the system can be described in the quantum jump 
approach through the no-jump Hamiltonian H — Hs + Hi. 
Here T-Ls consists of the decays and the free energy terms ll43l 
while. 


y|n)(ei| +iQi\ei){gi\b 
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H.c. 


( 1 ) 


Here the coupling strengths of the T** atom with the driving 
held and the incoming single photon pulse is respectively 
and Qi, while Vki is the van der Waals interaction among 
the Rydberg excitations of atoms k and 1. We solve the 
Schrodinger equation for the scattering stage assuming con¬ 
stant in Fourier space to hnd the rehection co-efficient. 
The (amplitude) rehection coefficient with the stored Rydberg 
excitation in atom k is given by TZkiuj) = (2KSfc(a;) — l) 



FIG. 1: Schematic outline of the phase gate (a) An input single 
photon pulse along with a driving Held induces a two-photon transi¬ 
tion to the Rydberg state |r) which is subsequently transferred to an¬ 
other Rydberg state \r'). Due to Rydberg interactions Vki among the 
atoms, other Rydberg states |r) within the range of the interaction 
potential, given by, the blockade radius of Vb, become off-resonant 
allowing no further excitation, (b) When an initial photon pulse is 
stored in the Rydberg ensemble, the second incoming photon cannot 
enter the cavity and is scattered off, which ideally induces a phase 
flip of TT on the scattered photons, (c) Dual-rail implementation of a 
CP gate (dotted box). A Bell state measurement can be implemented 
by combining the CP gate with Hadamard gates (HG). 


where. 
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Sk{uj) = \ K — iuj + 
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where the detunings are A; = We, — ojc, A/ — Ai + oj and 
6i = (wr, — We)—Wc with Iia;e,(r,) andrei(r/) being the energy 
and width of the excited (Rydberg) state |e)(|r)) in atom I 
and K is the cavity held decay rate. The rehection coefficient 
TZg{ijj) ll43l for no stored excitation is evaluated by setting 
Vki = 0 . 

To get an understanding of the scattering we study the be¬ 
havior of the rehection coefficient for resonant interactions 
6i = Ai = 0 and long lived Rydberg excitations (T^/ = 
0). For simplicity, we assume equal couplings and driving 
strengths on all atoms Qi — Q and Vti — VI (for the general 
case see mi Furthermore, if the photon pulse has a suitably 
long duration we can put w « 0 (see below). With these as¬ 


sumptions, we hnd from Eq. (S14i that TZg = 1 for no stored 


excitation; this is the perfect FIT condition. When an excita¬ 
tion is stored, the rehection coefficient becomes 


Ti-k = 


- 1 


where C* = Ct 
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tivity of the blockade, while C = \Q\‘^/{KVe) is the single 
atom cooperativity. Each atom I in the volume blocked by the 
|r^.) excitation, i.e. such that Vki ^ \Vl/ 2\^/Ve, will con¬ 
tribute with C in Cfc. On the other hand, those atoms for which 
Vki ^ 10/2 p/Eg will have a negligible contribution to Cb, 
and hence Cb gives the effective cooperativity of the blockade 
ensemble. The imaginary part depends on the shape of the 
interaction but for a uniform 1/r® interaction at resonance in 
a uniform cloud we find that \C'^\ = Cb ll43l (r is the distance 
between atoms k and 1). 

We now discuss the key feature of our work - the imple¬ 
mentation of a photonic CP gate via scattering from a Ry¬ 
dberg ensemble in either a single-rail or dual-rail encoding. 
The single-rail implementation uses the encoding discussed 
in the introduction and is shown schematically in Fig. [^(b). 
A first qubit is encoded in the vacuum and single photon state, 
|0) and |1) = J (ia;(/(w)ajje“*“*|0), respectively, of a first 
incoming pulse. This pulse is stored in the Rydberg ensem¬ 
ble such that the logical states |0) and |1) are mapped onto 
the ensemble being in the joint ground state |0) = |( 7 ^)| 0 ) 
and a Rydberg polariton |1) = ,r'i.)\0). This is 

achieved using the well established techniques of storage in 
atomic ensembles, which is known to have an error 1 /NC for 
any slowly varying pulse shape provided a temporally varying 
control field is used during storage EHH, followed by mi¬ 
crowave TT-pulse between jr^) and |r/). A second incoming 
photon pulse is then reflected from the cavity. This reflection 
can be from either an ensemble in the EIT configuration (en¬ 
semble in |0)), or from a blocked ensemble (|1)). As can be 
seen from Eq. ( |S17[ ) there is exactly a tt phase shift between 
the two situations: TZg = 1 for C/ = 0 and TZk = — 1 for 
\C* I ^ 1. Finally, the first stored pulse is retrieved from the 
ensemble. 

To evaluate the performance we calculate the Choi- 
Jamiolkowski fidelity of the gate. Since, in general, we have 
NC > Cb, the fidelity of the operation will mainly be lim¬ 
ited by the gate and we shall ignore imperfections during the 
storage. The fidelity can then be determined by ll43l , 


Erj ~ — 

CJ Jg 


duj\(l){uj)\‘^TZg{uj) 




(4) 


To account for errors due to imperfect Rydberg blockade, we 
evaluate the above fidelity and find 
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Here the third and fourth term are gate errors due to the fi¬ 
nite frequency width Aw^ of the incoming pulse. These terms 
arise predominately from the EIT bandwidth, which is much 
narrower than the variations of the blocked reflection coeffi¬ 
cient. For a narrow pulse Aw —> 0, the fidelity is only limited 


by the cooperativity of the blocked ensemble 1 — E/j cx 1/Cf,. 
Hence, as discussed in the introduction, it is the cavity en¬ 
hanced blockaded cooperativity, which is the main figure of 
merit for the gate. 

In the dual-rail encoding, both logical states |0) and |1) 
are represented by photons, but in two different paths. A 
schematic of the dual-rail CP gate is shown in Fig. [TJc). The 
first photon pulse in the upper two arms of the figure is first 
stored in a memory consisting of a Rydberg ensemble placed 
in each arm (for a polarization encoding such two memories 
might be realized by two different internal states of the same 
ensemble). A second photon pulse is then scattered from the 
Rydberg ensemble if it is in state |1) (upper rail in the figure). 
This scattering ideally induces a phase change of tt if there 
was a photon stored in the Rydberg ensemble, i.e. if both 
qubits were in state |1). As opposed to the single-rail imple¬ 
mentation, the dual-rail implementation has the possibility of 
conditioning on getting two photons in the output. Since the 
dominant error in the single-rail implementation is the loss of 
photons, this possibility allows for a substantial increase in 
the fidelity with only a minor failure probability of the gate. 
In view of a possible application of the gate for quantum re¬ 
peaters, discussed below, we consider the conditional fidelity 
of an EPR pair resulting from an entanglement swap realized 
with the gate using the full circuit in Fig. Neglecting 

again the error due to finite storage efficiency, we And that 
this fidelity is ll43l 


p — 
^ swap — 


/ dw|(/(w)p |2 -|-7^g(w) - J2k |cfcp7^fc(w)|^ 

16 Esuc 


( 6 ) 


where the success probability of the process is Psuc = 
/ tiw|(/(w)|2(2 + |E,(w)| 2 + |Efc|cfc|27^fe(w)|2)/4. Note 
that compared to Fq] in Eq. the only difference is due 

to the conditioning with a success probability Psuc < 1 and 
the way the mode function is treated. The latter is related to 
the fact that Eq. ( S19| l is the fidelity with a specific mode 
function. Keeping only the leading order contribution to the 
dispersion, we And the fidelity and success probability of the 
CP gate 
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Here the fourth term is again the leading order error from the 
spectral width of the pulse. In the limit of a narrow pulse 
Aw —> 0, we see that the conditional gate error 1 — Eiwap oc 
l/(Cb)^ for Cb ^ 1 is much smaller than for the single-rail. 
This comes at only a minor cost in the failure probability 1 — 
Psuc oc l/Cfe. The resulting dual-rail fidelities are plotted in 
Fig. I^as a function of the parameter Cb- For Cb « 8, the 
(post-selected) fidelity is found to be larger than 0.99 

In order to get realistic predictions, we use the experi¬ 
mental conditions of Ref. Il26l . with Fg « (27r)3MHz and 
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K « (27r)10MHz (corresponding to a finesse « 120) but 
a smaller beam waist wq = 15p,m. This gives a single atom 
cooperativity C = 0.025 and we take NC = 20 correspond¬ 
ing to a combined storage and retrieval efficiency of 90%. We 
assume a Rydberg line width 7 ^. = (27r)60 kHz ll44l corre¬ 
sponding to a coherence time of = l/ 7 r = 2.65 /rs (note 
that if the two ensembles in the dual-rail encoding are read out 
with the same laser the scheme becomes insensitive to phase 
fluctuations). With a pulse duration of T = l/Aw = 300 
ns and a driving strength of 11 = (27r)36 MHz the error due 
to finite bandwidth in Eq. ( |S22| i is below 2%. Taking the in¬ 
teraction V = (27r)8.31 • 10°/' MHz /rm® corresponding to 
two atoms with a Rydberg quantum number Ur = 90 and an 
atomic density of n = 0.25 one has Ct ~ 8.1 l43l 

which is sufficient to obtain high fidelities as show in Fig. 
Here, we ignore any effect of sample inhomogeneities, but this 
can be taken into account by suitable redefinitions of Cb and 

c'b m. 

F 



FIG. 2: Choi-Jamiolkowski fidelity (thin line) and post-selected 
swap fidelity (thick line) as functions of the blockaded cooperativ¬ 
ity Cb for a spectrally narrow pulse Alj —> 0. We assume |C/| = Cb- 

As a particular application of the gate, we consider long 
distance quantum cryptography based on quantum repeaters. 
We considered the ensemble based quantum repeater protocol 
proposed in Ref. 1461 , but replace the entanglement swapping 
with the procedure shown in Fig. [TJc). We calculate the secret 
key rate per repeater station as described in Ref. ia (assum¬ 
ing the distributed states to be Werner states) and compare the 
results to the original protocol (see Fig.|a. At the lowest level 
of the protocol, single excitations are stored in atomic ensem¬ 
bles using a Raman scheme and we include double excitation 
errors to lowest order similar to Ref. l46l . The performance of 
the protocol depends strongly on the repetition rate of this op¬ 
eration. Regardless of the repetition rate, the CP gate enables 
significantly higher communication rates since it allows near 
perfect Bell state measurements (for Ct 3> 1) whereas swap¬ 
ping operations based on linear optics have a maximal suc¬ 
cess probability of 50%. In Fig.[^ we also show the rate ob¬ 
tainable if single excitations are initially created perfectly and 
deterministically in the atomic ensembles, e.g., by exploiting 
Rydberg blockade 1291 . We find that for such a protocol, a 
cooperativity of Cb ^ 25 is sufficient to obtain a secret key 
rate of 1.5 Hz over 1000 km using 33 repeater stations. 

In conclusion, we have proposed an efficient method to im- 



FIG. 3: Secret key rate per repeater station (rsecret) as a function of the 
blockaded cooperativity (Cb) for a communication distance of 1000 
km. We compare the protocol of Ref. 1461 (Original) with a mod¬ 
ified protocol where the entanglement swapping is performed with 
the Rydberg CP gate (Rydberg). We consider an optimistic source 
repetition rate of 100 MHz and a more modest one of 1 MHz, as well 
as a perfect single excitation state created in the atomic ensembles, 
e.g., using Rydberg blockade 1291 . We assume an attenuation length 
of 22 km in the fibers and an optical signal speed of 2 • 10® km/s. The 
ensemble readout efficiency and photodetector efficiency are both as¬ 
sumed to be 90%. The steps in the curves reflect where the fidelity 
of the CP gate allows additional swap levels to be employed. 

plement a CP gate for photons. The gate combines the ad¬ 
vantages of cavity defined optical modes and cavity enhanced 
light matter interactions with the strong Rydberg blockade ob¬ 
tainable in atomic ensembles. As a direct application, the pro¬ 
posed gate can be used to improve the communication rate of 
quantum repeaters, but more generally the gate may serve as 
a building block for photonic quantum networks. 
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Supplementary Material 


THE REFLECTION COEFFICIENT 

The dynamics of the Rydberg ensemble in the cavity can be described through the no-jump Hamiltonian % consisting of the free 
energy and decay terms "Hs along with the interaction part as T-Lj = "HL-int + ^Ryd-int, where 


■Hs = ^ S(Ai - ire;)|ei)(e;| + h{5i - iTri)\ri){ri \ - ihnb^b 

i 

■HL-int = “X! - i^Wi\ei) {gi\b + }\.C. 

i i 

■HRyd-int = ^K\’ki\r'k){r'k \ ® k/)(n|, 


(S9) 


where the detunings Ai,6i, the linewidths Tg, T^., and the coupling strengths H;, Qi are as defined in the main text, while 2k is 
the cavity intensity decay rate. Note that all energies are measured relative to the cavity resonance, and hence the cavity term in 
Tis only involve the loss rate k. are as defined in the main text, while Vu is the van der Waals interaction among the Rydberg 
excitations of atoms k and 1. The incoming and outgoing photons are going to be accounted for by the input-output relations. 
After the storage of the first pulse, the wave-function of the combined field and ensemble with the initial excitation stored in |rj.) 
and one incoming photon is given by, 

1 ^') = XI / duJl3k{uj)ale-^‘^*\g^-^,r'k,0) + '^ i'^Ctk\g^~^,r'k,0) 

k k 

+ (SIO) 

i 

Here |0) is the vacuum state, Ceki and Crki are respectively the amplitude of being in the excited state |e) and the Rydberg 
state |r) when there is one stored Rydberg excitation in the /c* atom, while Cbk is the amplitude of the cavity excited state. We 
next evaluate the Schrodinger equation for the wave-function ( |S10| l together with the input-output relations to find the dynamical 
behavior of the amplitudes Cak, Ceki, Crki 


Cbk — ^ ^ CeklG* — KCbk + > 

I 

fl* 

Ceki ~ i^el)Cekl i^^Cj^kl GlCbki 

Crki — iA'rkl)Crkl “t” i^~^Cekl i^Crkl^kl- 

The outgoing field amplitude is then given by, 

/irM = 

where Cbk is found by solving the set of Eqns. ( |S1 ipO) using Fourier transformation. We thereby get, 

Cbk{u^) = V^/3L"(u;)Sfc(cc) 

\Gi? 


(511) 

(512) 

(513) 

(514) 


-1 


Sk{uj) = K — iuj + X! 


I + Tri+ksi+Vki-uj) 


(S15) 


Substituting Cbk into Eq. (S14i we get. 


TZk{i^) = 2k k — iuj + 


\Gi? 


I (Ee/ lAi loj) + 


- 1 . 


(S16) 


Assuming all fields to be resonant i.e. for <5; = A; = 0, a long lived Rydberg state (Tri = 0), and slowly varying photon pulses 
(w = 0) we get. 




j- i-i\ni/2\yvkire 


- 1. 


(S17) 
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which under the assumption of equal coupling strengths Qi = Q and Rabi frequencies fl; = fl, for the defined single atom 
co-operativity C — l^p/zcFe becomes, TZk = 2 (1 +C*)~^ — 1 , where C* = (l ~ *|f2/2p/Vfc;re). To get a simple 

physical understanding of the scattering dynamics, we shall first assume that all atoms are identical (homogeneous). We will 
consider what happens for an inhomogeneous ensemble in a later section. For the van der Waals interaction potential Vu = 
where r is the relative distance between the and atoms, we can evaluate C* with the sum converted to a 
volume integral —> / ndV. Thus we get for a homogeneous ensemble with an isotropic potential. 


C* = AnnC / dr j (1 + iCr°); C = 


(S18) 


JO 

We can write this integral as C* = Cb — iC^ and solved it to get, \Cb\ = \Cl\ = |(Cn7r^/'\/2C). Above, we have solved the 
scattering dynamics in the case where there was already a Rydberg excitation stored. In principle, we should also solve the 
dynamics without the first stored excitation. In this case, however, the excitations are completely independent of each other. We 


can then conveniently obtain the results for this situation by simple setting Vti = 0. Then from Eq. (S17i we get C* = 0 for a 
long photon pulse and hence TZg = 1. 

To investigate the effect of pulses of a finite duration, we now consider the bandwidth of the scattering coefficient. To do this, 
we perform a Taylor series expansion of the reflection coefficient about some central frequency ujq, 

1, 


TZkiuj) = 7^fc(a;o) + - Wo) +- Wo) . 

Here we have kept upto the second order in the expansion. The above three terms in the expansion are described by, 

2 . \ . Wo 1 


^/c(wo) = 


1+c; 


-1 


2i 


^ (1+Q)^’ 


(S19) 


(S20) 


du,Tlk\u:o = 




—4i 



with the parameters defined by, 

c: 


V 

1 

I ^ 


C 






1 + |n,/2|V(VH - uo)" 


cj. = cY. 


Fe Wo)Fe 

1 + |0;/2p/(Vfci — Wo)^ 


|SM! 


; 1 _ tJ_ i'“ti 

‘ i{Vk,-uJo)r, 

|Oj/2|2re/(Vfe;-Wo)3 


[■ 


|n,/2p 


' '<'rli{Vki-ojo)r, 


3 ’ 




[■ 


|Oi/2F 


' U + i(v„-wo)r. 


F. = cj;: 


|^^//2pwo /{Vki — Wo)" 




3 • 


(S21) 


(S22) 


(S23) 


(S24) 


(S25) 


Assuming the central frequency of the incoming pulse to be on resonance, we set wo = 0 and hence Eqs. ( S20[|S22 1 become 
substantially simpler and are described by. 


T^k = 

|a;o—0 — 

d'^T^k\ujo=0 = 


- 1 


1+C.. 

{1+C*) 


2 ’ 




(+% )' 

(1+Q)' 


r, r. 




(1+Q)® (1+Q)' 


— 4i 




(1 + Q)" 


(526) 

(527) 

(528) 
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Note that for the case of no stored photon in the ensemble, we have TZk —> Ti-g and one can also get the Taylor series expansion 
of TZg by setting Vki = 0. The expressions for such an expansion at resonance is the same as given by Eqs. (S26 - S281 but now 
with the set of parameters, C* —)■ C *—)■ C%,C*^ —>■ C*,C*^ —)■ C*, where the new parameters correspond to 
Eqs. ( |S23 - 1525 1 with Vki = 0. Erom the set of Eqs. ( S26|S28] l, we see that the leading order dispersive contributions are scaled 
down by a factor of (C*)^ and (C*)^ when we compare the situation with and without stored excitation in the first pulse. Hence 
the spectrally narrowest feature is the width of the EIT resonance without a stored excitation, and this will thus be the limiting 
factor for the bandwidth. 

We next analyze the behaviour of the parameters listed in Eqs. ( |S23|S25] ) in different limits of operation. We can find the 
blockaded part by considering the limits Vki ^ |^z/2p/re, while the contribution from the remaining EIT medium is found in 
the limit Vki ^ \^i/‘^\^/^e- To get a feeling for the expression in Eqs. ( S23|S25 i, we separate them into contributions coming 
from the blockaded atoms and that from the rest of EIT medium. 


Cb + iC'b — ^ 


C 


[■ 


|n/2|-‘ 


+ i 


E 


c 


vlr. 


, [1 + 


|n/2|4 

Vfcir^ 


(S29) 



P^Ol 

— 

ATOL p ^ e 


II 

* 

(S30) 

c* 

= cr 

r4 

c ® 

p^ 

= 0, 

(S31) 


where we have assumed all the Rabi frequencies to be equal such that H and C*“ = X); C/{1 -\- |H/2p/iVfc/re)^, = 
C*’' = + \Vl/2\'^/iVkiVef’, which scale as the number of blocked atoms Cb while N^jrp scale as the number of 

remaining unblocked atoms ~ N. Similarly, for the case of no stored excitation, we get. 


C* 

Cl 


= 0 , 


C = -NC 


|H/2|2’ 


CS=0; 


= -iNC 


\n/2\^ 


c* = o. 


(532) 

(533) 


Note that since in this case there is no Rydberg excitation blockade, only the EIT medium contributes and all the terms arising 
due to blockade are zero. 


CHOI-JAMIOLKOWSKI FIDELITY 

The Choi-Jamiolkowski (CJ) fidelity is a measure of how close two given quantum mechanical processes are. The idea is to 
apply each process to a particular entangled state and then calculate the fidelity between the two output states. Specifically, we 
assume that the two processes are described by the superoperators U and V. The superoperator U represents the ideal process 
that we want to accomplish and is assumed to be unitary. Hence, its action on some density matrix p can be be written as 

u{p) = UpU\ (S34) 

where C/ is a unitary operator. The actual physical implementation is represented by the completely positive trace preserving 
superoperator V. In general, it admits a Kraus (operator-sum) decomposition 

V(p) = ^ VipVi^ (S35) 

I 

with Vi^Vi = I {I is the identity operator). If we separate out the “no jump” evolution with the effective non-Hermitian 
Hamiltonian % in Eq. ( |S35[ ), we can write 

V[p) = VpV^ + Y, KipKj, (S36) 

i 

where V = ex-p{—i'Ht{/Ii) with tf being the time it takes to accomplish the wanted operation. The operators Ki form the Kraus 
decomposition of the part of the evolution where at least one quantum jump occurs. 

To find the CJ fidelity, we consider the superoperators X and X ®V that are tensor products of the original ones with the 
identity superoperator X. We pick an orthonormal basis set {| j)} for the d-dimensional Hilbert space that U and V act on. Now 
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we can define the state |$) = J2j |j) |j)/v^ that is an element of the original Hilbert space tensored with a copy of itself. Note 
that |<i)) is a maximally entangled state of these two copies. We will only consider a two-qubit gate, so that d = 4 in the above. 
After applying I and I onto the density matrix |$) (<i)|, we get a pair of new states 

pu = = (/(g)C/)|$)($|(/(g)C/t), (S37) 

pv = [I(g)V]|4>)($| = (7®y)|$)($|(J(g)y^)+^(J(g)Xz)|$)($|(7(g)if/). (S38) 

I 

The CJ fidelity is defined to be the fidelity of these two states. Since pu is a pure state, we get 

Fa = F{pu.Pv) = mi®U^)pv{I®Um 

= |($|(7® C/^y)|$)|2 + y]|($|(7® C/^7^z)|$)r 

i 

STORAGE AND RETRIEVAL 

The full physical process to implement the controlled-phase gate consists of storage of one photon, scattering of the second 
one, and retrieval of the first. The theory of storage and retrieval with an ensemble in a cavity is well established. In suitable 
regimes these results show that we have a mapping between a single mode of the atomic ensemble and a specific incoming 
or outgoing optical mode, and all other modes will be uncoupled IT]. Hence, the process of storage is described by a single 
parameter, which is the storage efficiency of a single incoming photon to create a specific spin wave 

\S)=Y,ak\g^-\r'^). (S40) 

k 

After scattering of the second photon, this spin wave will get multiplied by the reflection coefficient of the second photon such 
that it becomes 


\Sn) =^oikTlk{uj)\g^ \rfc), 

k 

where uj is the frequency of the second photon. Note here that the scattering coefficient may depend on which atom the first 
photon was stored in since different atoms may experience different degrees of blockade. For the retrieval, the cavity maps the 
particular spin wave ( |S40[ ) to a specific temporal mode. Hence, the amplitude of the retrieved photon is given by the shape of 
that temporal mode multiplied by the overlap 


{S\Sn) 

k 

In general, the retrieved wavepacket will also need to multiplied by the square root of the overall storage and retrieval efficiency, 
but we neglect this in our analysis. 

For the fidelity calculations, one would need to calculate the overlap of the retrieved photon wavepackets corresponding to 
I S') and \Stz). However, by the discussion above, the overlap of the photon wavepackets will be equal to the overlap of the 
spin waves |S) and \Sti). Hence, in the calculations below, we will directly calculate the fidelities by projecting the spin waves 
instead of analysing the retrieval. 


FIDELITY IN THE SINGLE-RAIL ENCODING 

In the single-rail encoding the computational basis is 

100(f)) = |5^)|0), 

101(f)) = \g^) I ‘10), 

110(f)) = y]afe|p‘^"\r;)|0), (S41) 

k 

|ll(f)) = J duj(j){u})ale-^‘^*\0). 

k 



10 


Note that this basis is time dependent due to the free evolution phase exp(—Hence, we define the ideal operation U such 
that it includes this free evolution phase. Specifically, if we denote the computational basis states at the initial time f = 0 by 
omitting the time variable, i.e. \jj') = \jf{t = 0)) 1), then the ideal operation of the controlled-phase gate is given by 


C/|00) = |00(tf)), 17|01) = |01(ff)), C/|10) = |10(ff)), C/|ll) = -|ll(ff)). 

Using the computational basis ( |S41| l we can write 

1$) = i (|00)|00) + |01)|01) + |10)|10) + |11)|11)) . 


Inserting this specific form of |<i)) into ( |S39| l we obtain 

1 


Fq] — — 

CJ 


(00|f/^U|00) + (01|C/1'U|01) + (10|C/^U|10) + (ll|f/V|ll) 


1 


— y (00|C/^Ki|00) + (01|t/^ifz|01) + {10\U^Ki\10) + {n\U^Ki\U} 
16 ■‘^1 


(S42) 


(S43) 


For the operators Ki, we assume that {jj'{tf)\Ki\jj') = 0, where j,j' G {0,1}. Physically, this assumption means that if a 
quantum jump (incoherent decay) occurs, the given basis will switch to another state of the physical system (possibly even one 
of the other basis states) but can never be driven back to the original state. I.e. if a photon is lost, it will result in a vacuum output, 
and thus is does not give an overlap with the original state. Under this assumption, we only need to compute the dynamics due 
to the non-Hermitian Hamiltonian. The detailed calculation is presented in Sec. [] In essence, the result is that the dynamics of 
the operator V can be described by the scattering relations 

uioo) = y)i0), 

U|01) = \g^) J du;TZg{u;muj)ale-^-^<\0), 

=y afciy"\T-fe)i0), 

k 

V\l\) = j Q^fc7^fc(a;)(^i(a;)a^,e"“*f|y"^r^)|0). 

k 


Gathering all the formulas in this section, the CJ fidelity becomes 

1 


Fa — 


16 


c^u;| 0 (w)p 7 ^g(a;) - / dw\(j){u})\'^^\ak\'^nk{u}) 


(S45) 


FIDELITY IN THE DUAL-RAIL ENCODING 

In this section we calculate both the CJ fidelity and the entanglement swap fidelity for the dual-rail encoding and show how 
they relate to each other. The circuit diagram of the entanglement swap operation is shown in Fig. 

The entanglement swap operation consists of evolution of the initial state (which is unitary in the ideal case) and a subsequent 
measurement. The evolution can be decomposed into a controlled-phase gate and Hadamard gates. If the Hadamard gates are 
assumed to be ideal, then the CJ fidelity of the whole evolution is equal to the CJ fidelity of the controlled-phase gate. We are 
going to use this fact in relating the CJ fidelity to the entanglement swap fidelity. 

The abstract definition of the CJ fidelity does not make any reference to a particular basis. In this section, in addition to the 
computational basis, we will also use the Bell basis, since it is the natural choice for the entanglement swap operation. The Bell 
states are 

10™) = 10+) = ^ (|00) + 111)) , 

|0™) = |0+) = ^(|O1) + |1O)), 

= \r) = ^ (|oo) - 111)), 

|0“) = |0-) = ^(|O1)-|1O)). 


(S46) 
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FIG. S4: (a) The circuit diagram of the entanglement swap operation. The numbers at the left edge indicate the label of the subsystem (qubit). 
In the circuit, the Hadamard gates are denoted by H and the controlled-phase gate is denoted by ip. 


In addition to the conventional names, we also give numbers to the Bell states, which will allow us to express summations in a 
simple way below. 


S4 


For the entanglement swap circuit of Fig. 

Bell pair |^^)34 between subsystems 3 and 4. Note that this initial state can be 


the initial state is one Bell pair 


l</>+>12 between subsystems 1 and 2 and another 
written as 


|0+)l2|0+)34 = J E I</>^'^'')i 4|0^'^")23 = |$)l423. 
1 J '=0 


This is exactly the state that is used as the input for the calculation of the CJ hdelity expressed in the Bell basis. After evolution 
of subsystems 2 and 3 as shown by the circuit and a measurement (the two detectors to the right), a Bell pair between subsystems 
1 and 4 is established. 

The practical implementation of the above circuit is shown in Fig. 1(c) of the main text. Whereas Fig. [^displays the extended 
four-qubit Hilbert space required for the calculation of the CJ and entanglement swap hdelity, Fig. 1(c) only displays the two 
central subsystems (2 and 3), but each of the two subsystems are represented by the photon being in two distinct modes. 

We dehne aj ^ to be the creation operator for subsystem 3 in state |0) with frequency w, and a| ^ to be the creation operator 
for state |1). For notational convenience, we dehne the states |0^)3 = Oq J\0) and |lcj )3 = d\ ^|0). Then in the dual-rail 
encoding, the computational basis is 


\0)2='^ak\g^ 

k 

| 1)2 = 

k 

|0(f))3 = J dw0(a;)e"*‘^*|O;^)3, 
|l(f))3 = J du}cj){u})e~"‘^*\luj)3. 


(S47) 


Here, \g^~^ ,r'i,)o are the states of the memory (the ensemble which does not interact with the second photon), and \g^~^ 
are the states of the cavity from which the second photon is scattered (see Fig. 1(c) of the main text). These two states correspond 
to subsystem 2 of Fig. S4 Subsystem 3 is encoded in photonic states which are not stored but only scattered. Note that all of the 
resulting computational basis states | 00 ) 23 , 101 ) 23 , 110)23 and 111)23 physically correspond to having two excitations. Hence, it 
allows for simple means of error detection: if less than two excitations are present at the end of the evolution, we know that an 
error has occured. In the dual-rail basis, the action of the operator V 23 that corresponds to the physical implementation of the 
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controlled-phase gate can be written 


V23 00)23 = E 

k 


^^23 01)23 = '^o:k\g^~^,rk)o j 

k 

[ c^w7^g(u;)())(a;)e"*‘^‘f|l^), 

1^23 10)23 = '^oik\g^~^J 

k 

[ du}4>{oj)e~"‘^'^‘\0k3, 

V23|ll)23= f duj'^ak'R-k{uj)(l){uj)e \rfe)i|l„) 

k 


(S48) 


For the operators corresponding to the full evolution of the circuit of Fig. we also need to describe the Hadamard operators. 
In the dual-rail encoding, the Hadamard operations are obtained by impinging the photons on beamsplitters which work on all 
frequency components separately. This is important for subsystem 3 (the scattered photon), since the frequency components 
will be multiplied with, in general, different reflection coefficients TZg{uj) and TZk{oj) depending on the input state. Hence 
the definition of the Hadamard operator here needs to be per frequency component, i.e. H3\0uj)3 = (|0tj)3 + |l^)3)/-\/2 and 
H3\1uj)3 = (|0(.j)3 — |la;)3)/'\/2. For subsystem 2, the single mode retrieval precludes any such difference in the mode shape 
for photons that are incident on the beamsplitters. Hence we can define the Hadamard operators to act on the spin wave states 
directly, H 2 \ 0 h = (| 0)2 + | 1 ) 2 )/V 2 and H 2 \l )2 = (| 0)2 - \ 1 } 2 )/V 2 .^ 

In analogy with Eqs. ( |S34| l and ( |S36| l, we define the superoperators IJ and V for the ideal and the real version of the circuit of 
Fig. S4 They can be written as 


^23 (P) — U 23 PUI 3 , 


and 

V23(p) = V 23 PV 23 + ^ ki^23pKj23, 

I 

where 

U 23 = {H 2 <8 H3)U23{l2 ® H 3 ), 

V 23 = {H 2 0 H3)V23{l2 ® H 3 ), (S49) 

ki ,23 — (H 2 0 H3)Ki^23{l2 ® H 3 ). 

Note that with the definitions ( |S42[ ), ( |S46| ) and ( |S49| ), it holds that ij23W^') = 

In this setting, not only the input states used for the entanglement swapping match the ones used for the CJ fidelity, also the 
actual operation itself has the same form: an identity operation acting on subsystems 1 and 4, while subsystems 2 and 3 are 
evolved according to either or V. The two output states are then 

= [X14 0W23] (|$)(^|) = {Iu®U23mmhi®Ul3)^ 

pp = [Ii4 0 V23] (mm) = (/l4 0 1^23)|^)($|(/l4 ® ^2^3) + 11(^14 ® ifi.23)|$)($|(/l4 0 i^/,23)- 

I 

Now we want to use the error detection property of the dual-rail encoding. We define the projection operators 


Pjj' =Ili® \j){j\2® 


MjDiL 


(S50) 


where j, / G {0,1}, and we also define their sum 


E 4 '- 

Jj'=0 


(S51) 


The projection operators of Eq. (|S50|l correspond to measuring the states |j /)23 on the detectors of circuit of Eig. S4 Note that 


for subsystem 3, we project onto the entire subspace that is spanned by the states |j^) instead of choosing a particular mode. This 
is equivalent to the assumption that all frequency components contribute to the probability of a “click” on the detector. On the 
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other hand, the operator of Eq. ( |S51| ) has a less clear physical interpretation. Formally, it projects a given state onto the subspace 
with two excitations. The motivation for defining such an operator is to be able to relate the CJ fidelity to the entanglement swap 
fidelity as we will see below. Since the entanglement swap fidelity can only be understood as a conditional fidelity (conditioned 
on the measurement outcomes corresponding to the operators of Eq. ( |S50| l), the CJ fidelity also needs to be conditional. 

Let us begin with the calculation of the entanglement swap fidelity. Using the states after the measurement has taken place 

^ ^ Pji'PyPjj' 

tr(4.pp4,)’ 

we can define the conditional fidelities for the entanglement swap 

F,,. =(4'|tr23(p,y0l4')i4. (S52) 

Here, we take the trace over subsystems 2 and 3, since the relevant question is how close subsystems 1 and 4 are to a particular 
Bell pair. The trace can be written as 

tr23(Pij')= y2 i tia;(n|2«|3Pij>L)3l»)2 = ■ [ duj{j\2{fj3Pv\fj3\j)2- 


Projecting pp onto \(jp^ )i 4 we obtain 


(4'|PvI<^''')i4 = ^4|4')23(4'|2343 + J E ^'^23 |4')23 (4' 

I 

For the dual-rail encoding, we have a stronger assumption about the operators Ki than for the single-rail encoding. We are going 
to assume that {n\ 2 {n'^\ 3 Ki\jj') = 0, where G {0,1}. Physically, this assumption means that the decay processes 

take the state out of the computational basis entirely, since any such decay will reduce the number of the total excitations to less 
than two. Then the expression for the fidelity (|S52|l becomes 


F ■/ = 


1 


4tr(4vpp4,) 


dui 


{j\2{jlj\3{H2 C) H3)V23{l2 C) Tf3)|4 )23 


with the trace in the denominator given by 

-I 4 /. 

tr(4'Pv4') = 4 E {n\2{n'j3iH2 0 P3)^23(/2 0 P3)|4')23 

n,n'—0 


For all j and / we get 


tr(PjyppPF,) = Yjj / duj\(l){uj)\^ I 2 + \TZg{uj)\‘^ + 


E |afcpPfc(w) 


and 


Pswap — Fjj' — 


16P, 


duj\cj){uj)\" 


2 + TZg(uj) — E! Wk\^'Ti,k{io) 


where we have defined the success 


probability = Ej,i'=o 


p 

/ 



j da;|(/)(w)p 

2 -f Pg(a;) p + 

1 

E |afcpP.fc(w) 

k 

) 


Now we look at the conditional CJ fidelity. Using the state 




tr(PppPt) 


(S53) 


(S54) 


(S55) 
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we can define the conditional CJ fidelity as = F{pjj, p'^). By the cyclicity and linearity of the trace, we have 

1 

tr(Pjj/pyPj^-,) = Psuc- 

jJ'=o 


The projection operator P has no effect on the states |$), hence F{pjj,Pp^P^) = F{pjj,p^), and the analysis reduces to 
finding the unconditional CJ fidelity and dividing by the success probability Psuc- The final result is 


^CJ = 


16Rl 


2 + duj\(l){oj)\‘^TZg{uj) - duj\(j){u})\'^'^\ak\'^'R.k{i^) 


(S56) 


Comparing Eqs. ( |S53| ) and ( |S56| l we see that the only difference is the order of integration and taking the absolute value. Thus 
in general, we have Pq < Pswap- If the bandwidth of the second photon is narrow compared to the the frequency variations 
of TZg and TZk, then the two fidelity measures become equal. The reason for this similarity is that the two measures consider 
the same input, but they do not consider exactly the same output. For the CJ fidelity the question we ask is what is the output 
with a particular mode, which we for simplicity take to be the same as the input mode. Possibly the CJ fidelity can therefore be 
increased by considering a more appropriate output mode. For the swapping fidelity we on the other hand consider everything 
which is incident on the photodetectors regardless of the temporal mode and hence this fidelity is higher. 


THE GATE FIDELITIES FOR THE RYDBERG CONTROLLED-PHASE GATE 


In the single-rail case, we evaluate the CJ fidelity ( |S45| l. To evaluate the quantity inside the modulus square we expand it and 
use Fq. (S19 1 and the corresponding expansion for TZg to get. 


ATZ = TZg — ^2 = 


2Ct, + 2iC' 


(1 + Cb) + iC'f^ 


(S57) 


ATZ' = 




= 2i 


1 NCT^ 


|fJ/ 2 |^ 


- 2i 


/V“ cv 

ISt 


2iCu 


|n/2|2 ;(1+Cf,)2 (1+Cfc)2 


(S58) 


ATZ" = 


TZ';^-Y,\ak\^TZ'l 








K |fJ/ 2|2 
Nlr^VlC 1 


{l+Cbf |0/2|4 


|T!/2|4 (1+Cb)3 ^Cl 


.1 


TT .. cr 


1 


(1+C,)3 


- 4f 


r 2 (i + c,)3 


(S59) 


In deriving the above expressions, we have assumed that the ensemble is homogeneous and that the potential is isotropic. Hence 
we have dropped the index k from Vki- We can then do the sum over k and given that ak are normalized, we have = 1 

in the above expressions. 

A closer look at Fq. (|S45|l suggests a further simplification which gives us. 


Fa = (4 + + 4Re[AP] + 2Re[AP"](Au;)2 + Re[A7^A7^"*](Aw)' 


(S60) 


where we have assumed narrow bandwid th of the p ulse and defi ned the variance of the incoming pulse as (Aw)^ = 
/ duj\(p{uj )\'^{lu — wo)^. Substituting Eqs. (S57 - S59l in Eq. (S 6 O 1 and assuming that Cb,C^ » 1 and < C^, we 

get. 


Fci — 


1 - 


(1 + Cb) 


il + Cb)^+Cl^ 


4[(1+C,)2+c;2] 2|H/2|4 


(A.r + 


(i + c.)^ + C 


1 / I NCTe 

2 |H/ 2|2 


C&(l + Ch)+Cj,^\ 2 


(S61) 
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Considering only the leading order contribution to the fidelity, we get, 


^CJ 


(1 + Cb) 

(i+oF+C 


NCTl 

WJW 


{Acor 


(I NCTe 




(S62) 


For the dual-rail case, we calculate a conditional swap fidelity ( |S53| ) and the success probability ( |S54| l. We can write ( |S53| l as 
1 


inswap - 


(4 -P |A7^p -h 4Re[A7^] + \ATZ'\^{Aujf -h 2Re[A7^"](Aa;)2 -p Re[A7^A7^"*](Aa;)2) (S63) 


which under the assumption that Cf,, 1 and Cj < C^, becomes. 


16Ra 


16 


1 - 


(1 -l- Cb) 


(i + a)2 + c;2 


SAcrL. „ Cbii + Cb) + c'/ 


(1+C6)2+c;2 |F!/2|4 


(Aw)" 1 -f 


(i+a)2+c;" 


1 NCT, 


K |11/2|2 

1 NCTe 




(i+a)"+c;" 


K |ll/2p 


(Aw)" 


1 , N^j^CTe 


\n/2\y \k |F!/2|2 ; (i + c6)2 


r(Aw)" 


(S64) 


The success probability ( |S54| l is then 
1 


p — _ 
j- sue - . 


(2 + |7^g|" + |7^;|"(Aw)" + Re[7^g7^"*](Aw)" + |7^fe|" + |7^),|"(Aw)" + Re[7^fe7^^](Aw)"), 


which under the assumption that C;,, C( 1 becomes, 


P — 1 — 

^ sue - 4- 


h ^ 

Cb 


NCTl 


(Aw)" 


(1+C6)2+c; 2 |r!/2|4' 

Substituting Eq. ( |S66| l and ( |S61| l into Eq. ( |S53| l, we get the expression for the conditional swap fidelity, 

Cb 


F ~ 1 — 

■*- swap — A 


3 C/ 3 

r 1 NCTe 1 
1 

2 

4[Cb + C'b^ + 2Cb] Pb + C/" + 2Cb] 2 4 

At |H/2|2_ 



i-f 


c2 + c:2 + 2a 


(Aw)" 


Einally, keeping only the dominant contribution to the gate operation, we get the conditional swap fidelity, 

,'2 on l2 


1 


p 1 - 3Cg-C 3 

[c2+c;2] 4 [c2+c;2]2 4 


1 ^ ACE, 


(Aa 


(S65) 


(S66) 


(S67) 


(S68) 


INHOMOGENEOUS ENSEMBLE 


So far, we have considered only a homogeneous ensemble without decay of the Rydberg level. In this section we discuss 
the case for an inhomogeneous ensemble. Eor simplicity, we only consider Aw = 0. Here the scattering dynamics depends on 
where the excitation was stored in the ensemble. Erom the fidelity expressions Eq. ( |S45| l and Eq. ( |S53| ) we see that the essential 
parameter is \ak\'^TZk which for <5; = A; = 0 is given by 


k k 


2 

1 I _ IPiP/Krei _ 


- 1 


(S69) 


When the ensemble was homogeneous, we defined the blockaded co-operativity Cb and C(, such that TZk = ■ Analogous 

to this for an inhomogeneous ensemble, we can define the blockaded co-operativity through J^k = 1/(1 -|-iCj,“^) 
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where. 


C inh 
b 


c 


^inh 

b 


Re 


Im 




Cl 


I^k i+|a/2|V(rHrei+iVfc,r,,), 


- 1 , 


J2k + J2i^k 


Cl 


i+|n/2|V(rHr.i+iVfc,rei), 


(S70) 


Note that contrary to the homogeneous case the above defined effective co-operativity for inhomogeneous ensemble also includes 
the effect of Rydberg decoherence on the scattering process. T hus, t o stu dy the Fidelity of the phase gate for an inhomogeneous 
ensemble and in presence of decoherence, the results in Eqs. (S62i and (S68 i can be utilized but now with Cb replaced by 
and by 
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